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ABSTRACT 



Wing rock at high angle of attack is an oscillatory 
lateral-directional motion phenomenon known to exist in 
some of today's high performance tactical aircraft. The 
motion has been consistently characterized as a lightly 
damped Dutch-Roll oscillation attributable to asymmetric 
wing stall. However, evidence gathered from wind tunnel 
simulations and at least one British study indicate that 
aerodynamic nonlinearities may be the source of wing rock. 
Regardless of the actual cause of the phenomenon, a study 
of wing rock has positive ramifications with respect to 
gaining a clearer understanding of the aerodynamics associ- 
ated with high angle of attack flight. This report presents 
the results of an investigation of wing rock which centered 
on the premise that two distinct nonlinear aerodynamic mech- 
anisms (aerodynamic hysteresis and a cubic nonlinearity in 
yawing moment) not only can cause wing rock but may drive 
it to a limit cycle oscillation as well. 
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I. INTRODUCTION 



Wing rock at moderate to high angle of attack is a char- 
acteristic known to exist in some of today's high performance 
aircraft, particularly those possessing the characteristics 
of thin swept wings and a slender fuselage. The terminology 
"wing rock" has been used to characterize an oscillatory 
lateral-directional motion exhibited by some aircraft when 
operating near stall, or when flow separation is evident by 
airframe buffeting. Although the motion has been observed 
and documented extensively during flight tests, its exact 
cause has remained essentially unknown. In an attempt to 
explain wing rock, fleet aviators, test pilots and aerodynam- 
icists have consistently advanced the theory which has gener- 
ally concluded that the apparent random rolling moments 
manifested by alternating wing rock motion at or near stall 
are due to asymmetric spanwise loading, i.e., alternating 
stall condition along the wing span. It is not the intent 
of this thesis to dispute such arguments but rather to present 
new arguments based on the premise that aerodynamic nonlinear- 
ities may not only cause wing rock but may drive it to a limit 
cycle oscillation. 

A. PAST EXPERIENCE WITH WING ROCK 

The author has had experience as a test pilot at the Naval 
Air Test Center and has personally observed what could be 
essentially described as a wing rock limit cycle while 
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conducting high angle of attack flight tests with the Northrop 
T-38A airplane. The T-38A is control limited longitudinally 
to the extent that it can not be made to achieve the so-called 
"classic aerodynamic stall", i.e., complete flow separation 
along the wing span. As a result, the normal 1-g stall is 
defined when a maximum angle of attack is reached due to 
control-limiting, moderate- to - heavy airframe buffet and a wing 
rock motion are evident with an average amplitude excursion in 
bank angle of _+ 20 degrees. This lateral motion occurs with 
essentially no heading deviation and will continue indefinite- 
ly, albeit accompanied by a large sink rate, until the angle 
of attack is reduced. The buffeting and wing rock will 
immediately cease and a complete recovery can be effected 
with very little additional loss of altitude, when afterburner 
thrust is applied. 

The question of whether wing rock motion in the T-38 (F-5 
series probably behaves the same) is indeed a limit cycle 
oscillation caused by aerodynamic phenomena other than asym- 
metric spanwise loading has never really been seriously con- 
sidered from a flight test standpoint. Instead, most test 
pilots and engineers associated with experimental/developmental 
flight testing quickly assume that an asymmetric stall condi- 
tion is responsible, since the motion is typically accompanied 
by separation as evidenced by airframe buffeting. Most swept 
wing aircraft exhibit wing rock motion to one degree or another 
near stall, although few can actually be said to exhibit motion 
behaving as a limit cycle. Since most aircraft are not control 
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limited, the approach to stall will often continue past that 
point where an oscillation similar to that exhibited by the 
T-38 might occur. The Vought A- 7 and F-8 series are notable 
examples of aircraft that can quickly reach a point where 
directional stability is lost, resulting in an almost immedi- 
ate departure from controlled flight, if the pilot does not 
respond quickly enough to reduce the angle of attack. 

Informal liaison with NASA Langley revealed that free 
flight model testing of the T-38A had on at least one occasion 
confirmed the existence of a wing rock limit cycle. Although 
credibility of this report would have been enhanced had it 
been possible to conduct an analysis of lateral directional 
motion at high angle of attack on a current production air- 
craft such as the T-38 or F5, the necessary information needed 
to conduct such an analysis was not forthcoming owing to 
proprietary considerations. Thus, it was necessary to inves- 
tigate and illustrate the concept of a wing rock limit cycle 
using the textbook airplane, F-94, as presented in Blakelock 
[Ref. 1] . During the course of the research the stability 
derivatives of the production F-94 were modified as necessary 
to simulate desired test conditions. To the extent that the 
F-94 served as a vehicle with which to demonstrate the concept 
of wing rock as a limit cycle, it was a satisfactory substi- 
tute in lieu of a current production aircraft. 

B. A NEW THEORY PROPOSED 

Explaining wing rock as an aerodynamic phenomenon apart 
from the theory of asymmetric wing stall poses some difficulty 
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partly from the fact that the motion closely resembles a 
lightly damped Dutch-Roll oscillation. Evidence obtained 
from closely controlled wind tunnel experiments [Ref. 2] and 
observations made during flight tests seem to suggest the 
possibility of a limit cycle oscillation wherein the amplitude 
and period of the motion may result solely from nonlinear 
aerodynamics. This' contrasts with the response of a lightly 
damped Dutch-Roll oscillation where the amplitudes are a 
function of the initial conditions. 

Some difficulty arises when attempting to identify an 
aerodynamic nonlinearity that may contribute to wing rock 
motion in wind-tunnel or flight tests since most dynamic 
analysis techniques, including aircraft parameter identifica- 
tion, are oriented towards a linearized model of the aircraft. 
Presumably, a nonlinear model could be assumed for parameter 
identification purposes but a question of uniqueness would 
arise since several completely different plant models might 
produce equally valid time history matches with actual flight 
data. Therefore, a significant step at this time would be 
the recognition of wing rock motion as a limit cycle oscilla- 
tion induced by aerodynamic nonlinearities. 

Ross [Ref. 3] , has shown that a wing rock limit cycle can 
exist when a nonlinearity of cubic form is present in either 
yawing or rolling moment with respect to sideslip angle. 

Ross used the Krylof-Bogoliubof method (KB method) with the 
assumption of slowly varying parameters to analyze the mathe- 
matical model of the motion described by the test airplane. 
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A portion of the analysis presented herein explores further 
the work done by Ross with, the objective of obtaining a 
simplified form of the theory she proposed as well as ampli- 
fying upon it. Preceding this, however, will be an analysis 
and discussion designed to support the claim that a wing rock 
limit cycle may also be the result of a nonlinear mechanism 
referred to as aerodynamic hysteresis. 

C. MOTIVATION 

The study of wing rock and its causes has direct applica- 
tion to an analysis of stall and post-stall gyrations which 
have long been difficult to predict and control on many 
tactical aircraft. The Vought F-8 and A- 7 are notable exam- 
ples. Increased knowledge in this area would surely enhance 
the ability to predict aerodynamic behavior, thus enabling • 
an earlier detection of undesirable flying qualities, or 
design improvement or, at least, a minimization of mission 
limiting or annoying characteristics. 
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II . METHOD OF INVESTIGATION AND EQUIPMENT USED 



The investigation of wing rock involved the use of the 
EAI C Electronics Associates Incorporated) model 580 analog 
computer as well as the Naval Postgraduate School's Depart- 
ment of Aeronautics digital computing facilities, namely, 
the Hewlett Packard, HP-9830. Mathematical models represent- 
ing the pertinent equations of motion were programmed on, and 
results compared from, both machines to ensure validity of 
programming. The output from these computing machines con- 
sisted of time histories of the motion exhibited by the test 
aircraft (F-94) when released from some initial condition. 
Output data were recorded by a graphical plotter or strip 
chart recorder attached as a peripheral device to the main 
computing units. These time history plots were used to make 
both qualitative and quantitative assessements of the motion, 
including measurements of critical damping ratio, frequency 
and/or period. Digital as well as graphical output was ob- 
tained using the HP-9830. Time history plots appearing herein 
were initially made using the analog computer, but final copies 
were obtained from the HP-9830 plotter as it produced a more 
compatible form for this report. The analog computer was 
superior to the digital machine in terms of its ease of pro- 
gramming or modification of the program during execution where 
it was possible to view readily the effects of altering certain 
key aircraft parameters. 
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The mathematical model used for laboratory tests consisted 
of the linearized lateral-directional equations of motion with 
certain simplifications to be presented later. As mentioned 
earlier, the F-94 airplane as presented in detail in Blakelock 
[Ref. 1] was used as the test vehicle for conducting laboratory 
simulations of wing rock motion. In its production configura- 
tion the airplane exhibits no inherent instabilities with 
respect to the lateral-directional modes of motion. Modifica- 
tions were made to its production stability derivatives as 
necessary to conduct desired simulations. 
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III. RESULTS AND DISCUSSION 



The discussion will focus on the particular lateral- 
directional motion exhibited by some aircraft within their 
respective flight envelopes and referred to hereafter as wing 
rock. The first part of this section will be devoted to a 
presentation of the analysis and investigation of wing rock 
as a limit cycle induced by aerodynamic hysteresis. The 
second part will present a similar analysis and investigation 
of wing rock as a limit cycle induced by a natural instability 
in the aircraft that is constrained by a nonlinearity in one 
or more stability derivatives. The reader is advised to refer 
to pages 8 and 9 for definitions of symbology used throughout 
the following discussion. 

A. THE GENERAL LIMIT CYCLE OSCILLATION 

A limit cycle oscillation is similar physically to the 
motion described by the typical undamped harmonic oscillator 
represented by: 

X -f- ^ X = O 

which has solution: 

Xft) = Si S/a/ ( ) !Z< t t- 4>) 

X(t) - Cos (V^T i: -h 

The motion amplitude of the linear undamped system will depend 
upon the initial conditions, and hence different initial con- 
ditions will result in different amplitude histories. 
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The linear undamped system will, of course, exhibit no damping 
behavior of the motion. By contrast, a nonlinear system that 
demonstrates a limit cycle behavior will have the same ampli- 
tude description for motion regardless of the initial condi- 
tions. A system that is weakly nonlinear with a limit cycle 
behavior will appear to have the properties of a linear system; 
however, the characteristic of motion amplitude regardless of 
initial condition is a distinctly different attribute of the 
nonlinear system. One should recognize that the above differ- 
ences are somewhat heuristic and a mathematician would char- 
acterize the distinctions with considerably more detail. 

B. AERODYNAMIC HYSTERESIS 

Wind-tunnel tests on aircraft models often show that the 
rolling moments measured on a symmetric model, when yaw angle 
is zero and angle of attack is increased, are approximately 
zero-valued until the angle of attack approaches close to 
stall. Rolling moments tend to be erratic at or near stall, 
and it is customary to attribute these traits to slight 
asymmetries in the model or tunnel flow, with a resulting 
early stall of one wing relative to the other. A recent wind- 
tunnel test [Ref. 2] upon a sting-mounted model having a strain 
gage balance with appropriate electronic filtering, used the 
test technique of conducting a slow, continuous yaw angle 
sweep through a complete cycle of yaw motion. When the cyclic 
yaw sweep was conducted at or near stall angle, both rolling 
and yawing moments showed incremental shifts in value that 
appeared to be related to the direction of the slow yaw sweep. 
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The curves of rolling moment, when plotted versus yaw angle, 
gave evidence of aerodynamic hysteresis. 

Simply stated., aerodynamic hysteresis is a type of behavior 
wherein an aircraft's response to a given input is such that 
it does not follow the same path returning to equilibrium as 
it did when initially displaced. Of importance here is the 
fact that the aircraft remains positively stable. Figure 1 
is an illustration of stall hysteresis wherein the airplane 
experiences a sudden decrease in lift coefficient as a result 
of complete flow separation along the wing span at stall angle 
of attack. The airplane recovers to a lower lift coefficient 
and angle of attack only to repeat the cycle if the controls 
are held in a pro-stall condition. This type of motion is 
often referred to as "stall porpoising" and is common to most 
aircraft with high aspect ratio wings. This example of longi- 
tudinal hysteresis is presented to the reader as an aid in 
understanding the discussion to follow pertaining to lateral- 
directional hysteresis. 




Figure 1. Stall porpoising and longitudinal 
hysteres is . 
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1 . Theory 



The analysis presented here will include some simpli- 



fications where deemed appropriate. The interest is to use 
state variable format to represent the mathematical model of 
the pertinent lateral-directional motion. The state variable 
format will be extended in an illustrative sense to demonstrate 
a wing rock limit cycle due to aerodynamic hysteresis. It is 
assumed that the lateral -directional motion of the aircraft 
occurs in such a manner that the airplane remains on a straight 
flight path, hence yaw and sideslip angles are related by: 



Further, it is assumed that the aircraft is operating at a 
fixed angle of attack with all control surfaces fixed. 
Therefore , 




( 1 ) 



A u- z yawing moment 




(02a) 



and 




moment 




(02b) 



Since A = = - j3 and A = - equations (02a) 



and (02b) can be simplified to: 



/ = +K/? 



(03a) 



and 




(03b) 
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Further modifications have been made yielding a dimensional 
form of the linearized lateral-directional equations of motion 
The assumption that a straight flight path is maintained 
effectively eliminates any influence from sideforces and the 
usual accompanying spiral mode. Such an assumption is reason- 
able and consistent with flight test observations of wing rock 
motion. 

The freely vibrating system represented mathematically 
in equations (03) may be cast into state variable format as: 

H-HH 

which corresponds to the shorthand notation for the following: 




o / o 
-/> -N f 



- ft ^A ^ . 



‘/-J 




(04) 



A third equation or identity ft-/? is necessary to accommo- 
date all the state variables in the matrix format. 

To digress briefly, J x] represents the state vector 
which contains the components of sideslip angle, ft , sideslip 
velocity, ft, and roll rate,y^. The plant matrix [ ] is of 
order 3 and contains the dimensional stability derivatives 
(standard NASA convention) which resulted from the previous 
analysis. These derivatives are explicitly defined as follows 

= o 



n 



= * 



L l/ 



= -ftf = ^ 
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Equation (04) represents an eigenvalue problem of a freely 
vibrating system in which the eigenvectors are the modes of 
motion and the eigenvalues contain the corresponding char- 
acteristics of damping, natural frequency, and period. In 
its present linear form the system described by equation (04) 
has a solution consisting of three roots (one real and one 
complex conjugate pair) corresponding to the roll convergence 
(subsidence) and Dutch-Roll modes respectively. Again, making 
use of the state variable format the solution to the freely 
vibrating system may be expressed as: 



the vector of values for the state variables at time equal to 
zero. In essence, the state transition matrix is used for 
modal analysis to project or propagate the aircraft's motion 
numerically from some starting values (initial condition) to 
another set of values at a specified time greater than zero. 

The equations of motion for the dynamic system with 
control applied is given by: 




(05) 



where [^(^)] is the state transition matrix and 7 X ( 0 ) f is 




( 06 ) 
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where the control matrix [S] contains the control power 
derivatives, for yawing and rolling moments associated with 
the aileron and rudder controls. The vector fa J contains 
the magnitudes of the inputs for the respective controls. 

The general solution to equation (06) is expressed as: 

r t 

fm) - [*<#] fx®} + I fe(-T-)][e]{umJ jt- to?) 

where the integral term is a convolution integral. This is a 
general solution in the sense that it applies to any of the 
standard input methods, i.e., step, ramp and impulse. The 
analysis and experimentation to which this section of theory 
applies was restricted to the case of step inputs, which is 
probably the most common control input method used in aero- 
dynamics for studying aircraft dynamic stability. 

The specific solution to equation (06) for a step 
input of control is expressed as: [Ref. 4] 



|x(*)J - feco)J +fe] few - l][ d] fej 



(08) 






where [/?] and [I] are the plant inverse and identity matrices 
respectively . 

Figure 2 contains a representative form of the aero- 
dynamic hysteresis assumed in the analysis of rolling moment. 
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Figure 2. Assumed form of roll hysteresis. 

The motion of the system corresponds to a periodic variation 
of the state variables such that the amplitude of the sideslip 
angle is assumed to be bounded by: 

Referring to figure 2, the hysteresis loop is defined by the 
dashed lines. The value of the rolling moment corresponding 
to a particular sideslip angle will be the sum of the produc- 
tion aircraft's linear roll-due- to sideslip and the additional 
hysteresis term, which has constant magnitude but varies in 
sign according to the sign of the sideslip velocity. This 
latter condition alludes to the manner in which the increment 
of rolling moment was applied, i.e., through a relay action 
mechanism ( + ) driven by the sign of sideslip velocity. Side- 
slip velocity will naturally be plus or minus in sign according 
to the direction of yaw, hence the increment of rolling moment 
was added or subtracted accordingly. Thus, the rolling moment 
due to sideslip as shown in Fig. 2 may be expressed as: 

Roll due to sideslip = L (&) = ^ + 5/<£/V (/&) 
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Figure 2 should not be construed to represent the only form 
that roll hysteresis may take but was used in this analysis 
to illustrate the concept. Nothing has been stated explicitly 
or implied up to this point with respect to the direction 
assumed by the hysteresis loop as the aircraft oscillates in 
yaw, i.e., clockwise or counter-clockwise, but this will be 
clarified later in the presentation of experimental results. 
Although not shown, yaw moment hysteresis may be assumed to 
be of similar form as that presented in Fig. 2. 



was necessary to solve equation (08) after some simplifications 
required to facilitate the ease of computation were implemented. 
Under the assumption of a limit cycle, the motion is periodic 
and the amplitude remains constant, which can be expressed 
in part as : 



where T is the period. Referring back to the column vector 
of state variables in equation (05) and with a convenient 
choice of the time origin: 



This is a necessary and sufficient condition to assure proper 
establishment of the relay action on roll moment hysteresis. 
The forcing function portion of equation (06) may be expressed 
as : 



In order to predict the behavior of the aircraft it 




X z (0) = (T/z) = X(T) = /?= o 




( 09 ) 
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where 



AL 




The value of unity in vector element u(3) of equation (09) is 
used to denote a step input of the roll moment increment, AL . 
This can be equated to a step input of aileron by the pilot, 
though roll hysteresis should not be considered as having that 
as its source, in keeping with the requirement for this analysis 
that controls be -Fixed. An external source may, in fact, be 
the origin of roll hysteresis. Carrying the previous logic 
further, it is valid to consider the motion over a half-cycle 
since periodicity once again requires that: 



Substituting these simplifications into equation (08) yields: 



Equation (.09) can be solved for the period T by numerically 
iterating for the value of T/2 which will cause X2(0), side- 
slip velocity, to vanish. However, it is necessary first to 
solve for the state transition matrix [^(tb)] in order that 
all values are known except T/2. Coincidental to computing 




and 



XjW - - *i(TA) =o 
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the period, the maximum values of the state variables X^(0) 
and X-^Q) (sideslip angle and roll rate, respectively) were 
determined, since these necessarily occurred at the state 
point where X ? (0) vanished. Computationally, this was a 
simpler method than integrating equation (07). A sample 
solution will be presented later. Upon examining equation 
(09) it should be noted that the period of the limit cycle 
in no way depends upon the sign of roll moment hysteresis; 
however, the amplitude of the state variables is a linear 
function of the magnitude of roll moment hysteresis, AL- 
The same procedure as outlined above could be applied to an 
analysis of yaw moment hysteresis. 

2 . Analysis 

For reasons given earlier in the introduction, the 
F-94 airplane was chosen as the test vehicle for conducting 
laboratory simulations of wing rock motion. In addition to 
being in the landing configuration, the following pertinent 
data relative to the F-94 was obtained from Blakelock [Ref. 1] 

Cj A = .0278 RAD" 1 



U = 135 mph = 198 fps 
q = 46.6 psf 
S = 239 ft 2 
b = 37.3 ft 

lx = 7.16 x 10"^ slug-ft 2 



< 4 . - -- 0916 " 

C 



*10 



= .105 



Iz = 33.01 x 10 -5 slug-ft 2 

% 

Cjf = -.450 



= -.048 7 RAD" 1 



Cft jp = ", 0 5 3 
= - . 210 



For these values the plant matrix of equation (04) became: 
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Q 



1 



0 



[A] - 



-1.3214 .2491 .0629 

-2.822 -1.517 -2.4557 



( 10 ) 



a. Numerical Predictions 

The system of three first order, ordinary linear 
differential equations in (04) were solved by using the .com- 
puter program subroutine BASMAT, listed and described in Melsa 
and Jones [Ref. 5] . The Standard Fortran IV code was trans- 
lated into the BASIC program language code for adaptation of 
BASMAT to the HP-9830 digital mini-computer. The eigenvalues 
of equation (04) using the plant defined by equation (10) were: 

X 1 = a = -2.4475 SEC -1 
X 2 3 = b + ic = -.12873 + 1.1755iSEC" 1 



These values contain characteristic information pertaining to 
the two modes of motion which the aircraft exhibited as a 
freely vibrating system, i.e., no hysteresis imposed. These 
two modes were the roll convergence and Dutch-Roll modes. The 
roll convergence mode is indicated by the single real root, A^, 
while the Dutch-Roll mode is represented by the remaining 
complex roots, A 2 3 * This information was needed for compari- 
son with the motion solutions obtained when hysteresis was 
applied (to be presented later) . The Dutch-Roll mode indicated 
the following characteristics with respect to period, frequency 
and damping of the stable motion. 
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Hence : 



0.12873 sec 



- 1/7 17 5 



/V 



1.755 sec 



-1 



5 = 0.109 



LU 



1.183 sec 



-1 



n 



T 



d 



5.345 sec . 



These calculations show that the Dutch-Roll mode is indeed 
stable, as expected for the production aircraft, with a damp- 
ing level and undamped natural frequency of 10.9% of critical 
and 1.183 rad/sec, respectively. The roll convergence mode 
was not considered further, since the primary interest was in 
the oscillatory mode of motion. 

The state transition matrix, [0(f“)], was obtained 
likewise from BASMAT and is presented below. 




( 11 ) 



where 



e 0.0114 0.0012 0.0093 

[R] = -0.0279 -0.0030 -0.0228 

1.2146 0.1318 0.9916 



Q . 9886 -0.0012 -0.0093 



0.0279 1.0030 0.0228 



-1.2146 -0.1318 0.0084 



0.1320 0.8531 0.0184 



-1.1791 -0.1084 0.0086 



-0.0050 -1.0305 -0.0236 
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and a, b, c derived from = a 

X 2,3 = b - 1C 

With this information input to equation (09) , and applying an 
iterative solution technique, the period and amplitudes of the 
state variables for the cases of imposing roll and yaw moment 
hysteresis either independently or concurrently was obtained. 

The results for varying values of roll and yaw moment hysteresis 
are presented in Table I. 



TABLE I 

Results of Numerical Solution of Eq. (09) for 
Different Magnitudes of Roll § Yaw Hysteresis 



Type 

Hysteresis 


Magnitude 

of 

Hy steres is 


Max P 

(Radians) 


Max 

(Rad/ sec) 


Period T 
(sec) 


Roll 


1 


. 0922 


. 2948 


5.5977 


Roll 


2 


. 1844 


. 5897 


5. 5979 


Yaw 


.05 


. 2105 


- . 2551 


5.3459 


Yaw 


.1 


. 4210 


- . 5102 


5. 3457 


Roll 


-1 


- . 0922 


- . 2948 


5.5977 


Roll 


(1) 


1 


.3007 


. 0419 


5.4219 


+ 

Yaw 




.05 


Roll 


(1) 


. 05 


.1714 


-4. 514E-03 


5.4121 


+ 

Yaw 




. 5 



Cl) Only condition for which roll rate, -f? , is not a 

maximum. All other values are necessarily maximums. 
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The iterative solution technique of interval 
halving, as described by Gerald [Ref. 6] was used to solve 
Eq. (09) for the values in Table I. Two function values of 
opposite sign were initially obtained by inputing estimates 
of the period. Thus started, the interval halving routine 
converged in an iterative fashion to a solution. The computer 
program used for this purpose and implemented on the HP-9830 
is presented in Appendix A. 

The state transition matrix, with its linear com- 
bination of the eigenvectors (mode shapes) , allowed identifi- 
cation of the individual parameters comprising the Dutch-Roll 
mode. This identification follows from equation (11) and the 
following trigonometric identity: 

S/A/ 9 +6 CO 5 e = K (ir & *• -f cas ») = KM (e + f) 

where 

2 2 

K = A + B = magnitude of motion vector in polar 
coordinates 

= Arctan B/A = phase angle of motion vector in 

polar coordinates. 

Applying this mathematical form to the second term on the 
right of equation (11) results in: 

[fj.. C*S (Ci) + fc ] £>i Sin (Or) = « M <9 + -f- 







f 


\ 


, 1/2 


where 


K i,: 






j = magnitude of motion vector 


and 


'i.j 


= Arctan 


[F] iJ ' 


= phase angle of motion 






vector 



i,j = array locations. 
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Using this and appropriate initial conditions such that 





/ 



/ 

o 

0 



the resulting relative magnitudes and phase angles of the 
elements comprising the total Dutch-Roll motion vector when 
normalized with respect to sideslip angle were: 



I 




ft = 1.00 

/? = 1.183 
•f = 1.218 
<p ~ 1.030 




f/5 = 0.0 Deg. 

= 96.25 Deg. 

fV = 187.38 Deg. 

fp m 91.13 



Although bank angle was not a state variable, its 
value can be determined as a function of time according to the 
relationship fifty ~ y?(t) dt , or approximated by fifty ? • 

The integration approach was used on both the analog and digital 
computers. Bank angle was introduced to gain a better under- 
standing of the motion's character, specifically with respect 
to the relative amplitude ratio of bank angle to sideslip angle. 
A positive value for the phase angle above indicates phase lead, 
b. Verification of Numerical Predictions 

The numerical computations presented in the fore- 
going were verified qualitatively by modeling the subject air- 
craft in the laboratory on the EAI-580 analog computer. The 
circuit required for the hysteresis model is shown in Appendix 
E. The same model using the system of equations given in 
equations C.03) was programmed on the HP- 9830 for the purpose 
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of obtaining quantitative data. Digital propagation of air- 
craft motion was obtained by a numerical time-solution of the 
differential equations using the Euler Predictor- Corrector 
method given in Ref. [fc] . The computer program used on the 
HP-9830 for this purpose is presented in Appendix B. A 
separate program (Appendix C) for plotting output from the 
program in Appendix B was necessary to avoid overflowing 
available memory capacity due to large storage arrays used. 

Time history plots of aircraft response are pre- 
sented in figures 3 through 10. These computer generated 
plots were obtained by numerically integrating, as previously 
mentioned, the following governing equations of motion: 

/$ - - /S -h /v^ J - z 7 *■ An * jy<jvy (/?) 

~ J-A. AL Yr J7<f/V ) 

These are equations (03) with roll and yaw hysteresis terms 
added. The third equation contained in matrix equation (04) 
was merely an identity and thus served no useful purpose in 
performing numerical integration. 

The time history responses in figures 3 through 6 
are arranged to facilitate comparison of the characteristics 
of the two forms of motion, i.e., the damped Dutch-Roll and 
limit cycle oscillation. The Dutch-Roll mode was not dupli- 
cated on all plots; therefore, it will be necessary for the 
reader to refer back to figures 3 through 6 for comparison 
purposes. All of the time history data were plotted with 
sideslip angle as a reference parameter. 

Figures 3 through 6 show the limit cycle resulting 
from the application of an increment of roll moment hysteresis 
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RRD IRN5 UR RRD IRN5/5EC RRDIRN5 DR RRDIRN5/5EC 




TIME < 5EC ) 




TIME (SEC) 



Figure 3. Stable Dutch-Roll (Upper Plot); Wing Rock Limit 
Cycle Induced by Roll Hysteresis (Lower Plot) , 

3 and 3 Comparison. 
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ld oa 3- m 



TIME C5EO 




TIME C5EC) 



Figure 4. Stable Dutch-Roll (Upper Plot); Wing Rock Limit 
Cycle Induced by Roll Hysteresis (Lower Plot), 

$ and p Comparison. 
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rN 
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T* 



03 



time: cseo 




Figure 5. Stable Dutch-Roll (Upper Plot); Wing Rock Limit 
Cycle Induced by Roll Hysteresis (Lower Plot); 

$ and ^ Comparison. 
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Figure 6. Phase Plane Plots of Dutch-Roll and Roll 
Hysteresis Induced Wing Rock Limit Cycle. 
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(^\L= 1) via the relay action mechanism. The relay action was 
simulated on the analog computer by patching the output of 
sideslip velocity through a comparator and relay module whose 
output of +_ 10 volts was a function of the sign of /S’. The 
same mechanism was simulated on the digital computer by making 
use of the predefined function, SGN, which returns a plus or 
minus 1 value according to the sign of its argument. 

Two rather unusual features that require some 
elaboration are present in the time history plots. Specifi- 
cally, the unusual shape of the roll rate curve with roll 
hysteresis applied was believed due to the sudden application 
of a rolling moment step input via the relay action. The 
lower of the two peaks per half cycle was judged the more 
accurate for analysis (refer figure 4). Secondly, a solution 
was obtained to equation (08) using a negative value of the 
roll hysteresis, which except for a change in sign was identical 
to the solution obtained using a positive value of equal mag- 
nitude roll hysteresis (Table I). However, both the analog 
and digital computers gave a time history response for the 
same value of hysteresis that was heavily damped or essentially 
deadbeat after one cycle. The numerical prediction was there- 
fore assumed to be in error, although this did not reflect on 
the validity of the mathematical modeling used. Rather, it 
was concluded that the predicted solution failed to account 
for the overall interaction of the production aircraft's sta- 
bility derivatives during a dynamic propagation of the motion. 
Stated another way, the hysteresis action would only traverse 
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in the counter-clockwise direction about the loop in figure "X. 

A similar check was made on yaw hysteresis with identical 
results . 

Figure 7 shows the limit cycle caused by applying 
yaw hysteresis via the same relay mechanism used for roll 
hysteresis. When compared to the corresponding time history 
plot for the Dutch-Roll oscillation in figures 4 and 5, the 
motion induced by yaw hysteresis in figure 7 bears a remarkable 
resemblance. The phase angle relationship and period (Table I) 
are identical to the Dutch- Roll. From this, one may conclude 
that yaw hysteresis is essentially a forced Dutch-Roll response. 

The combination of roll and yaw hysteresis reflects 
to some extent the relative influence of each on the total 
aircraft response. As one might intuitively expect, the result- 
ing motion was a compromise or average of the individual effects 
of roll and yaw hysteresis as previously seen. In Table I it 

was pointed out that the maximum value for roll rate was not 

* 

necessarily achieved when sideslip velocity,^, was zero. This 
was confirmed by the plot in figure 9. (Although not explicitly 
shown in figure 9,(3 was obviously zero when / was an extremum). 
There were many more combinations of roll and yaw hysteresis 
that could have been tried; however, this was not considered 
feasible nor productive. Again, the relay action can be seen 
in the behavior of the roll rate response in figure 9. The 
distortion of the roll rate response curve was due to the 
combined effects of roll and yaw hysteresis. Likewise, the 
phase shift for bank angle is a compromise between that seen 
in figures 7 and 5. 
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Figure 7. Time History of Divergence to Wing Rock 
Limit Cycle Due to Yaw Hysteresis. 
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Figure 8. Phase Plane Plot of Divergence to Yaw 

Hysteresis Induced Wing Rock Limit Cycle. 
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Figure 9. Wing Rock Limit Cycle Due to Combined 
Roll and Yaw Hysteresis. 
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With the exception of figures 6 and 8, the limit 
cycle oscillations were not completely developed at the ter- 
mination of the time allotted for the numerical integration. 

This was due primarily to the limited plotting space avail- 
able, whereas in figures 6 and 8 the phase plane plot was not 
so restricted, and clearly shows convergence or divergence to 
the limit cycle amplitudes after release from the initial 
conditions . 

For the sake of visualization, one could equate 
the values of roll or yaw hysteresis to a deflection of the 
appropriate control surface, i.e., a nominal +_ 10% of avail- 
able control power from the aileron or rudder control surfaces. 
From the analysis conducted on the analog computer, with its 
real time ability to reflect changes in the problem variables, 
it became apparent after some experimentation that the cross 
coupling between roll rate and yawing moment through the 
adverse yaw stability derivative , N-f , though weak, was respon- 
sible for excitation of the limit cycle when roll hysteresis 
was applied. Yaw hysteresis, on the other hand, acted directly 
through the directional stability derivative,^, to produce 
a yawing moment. 

3 . Summary of Hysteresis Induced Wing Rock 

There are several significant points worth considering 
from this analysis and the data presented in Table I and 
figures 3-9. These include: 

Cl) The amplitudes of the limit cycle were directly dependent 
upon the magnitude of applied hysteresis. 

(2) The period of the limit cycle was independent of the 
magnitude of applied hysteresis. 
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(3) While the roll rate was 180 degrees out of phase with 
the sideslip angle for the Dutch-Roll oscillation, it 
was exactly in phase for the limit cycle due to roll 
hysteresis . 

(4) In addition to phase differences, the ratio of bank 
angle to sideslip angle was three times greater for 
the limit cycle due to roll hysteresis than for the 
Dutch-Roll oscillation. 

C5) The limit cycle for yaw hysteresis is essentially a 

forced Dutch-Roll response. Frequency and phase angle 
relationships are unchanged from the Dutch-Roll. 

(6) Reversing the relay action dependence on the sign of 
8 results in a deadbeat response. Therefore, a wing 
rock limit cycle is dependent upon the sign of roll 
or yaw hysteresis. 

(7) Combined yaw and roll hysteresis produces a limit 
cycle response which is essentially a compromise or 
average of the two effects. 



C. ANOTHER FORM OF AERODYNAMIC NONLINEARITY 
1 . Previous Research 

As stated earlier, some research has been done and 
reported by Ross [Ref. 3] on wing rock as a limit cycle caused 
by nonlinear aerodynamics interacting with an unstable Dutch- 
Roll oscillation. Ross reports that a wing rock limit cycle 
occurred on the Handley Page 115 testbed aircraft as the 
result of increasing the angle of attack to a critical point 
where any small lateral-directional perturbation excited an 
alternating wing rock motion, which steadily grew to a con- 
stant amplitude with the controls fixed. The motion would 
immediately damp to zero when the pilot commanded a reduction 
in the angle of attack below critical. The Handley Page 115 
is a British aircraft designed to serve as a flying testbed 
primarily for research on the Supersonic Transport and has 
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geometric characteristics similar to many current tactical 
aircraft, i.e,, thin swept delta wings and a slender fuselage. 

In modeling the motion observed in flight and wind- 
tunnel tests, Ross essentially took the linearized equations 
describing lateral -directional motion and imposed a nonlin- 
earity on yawing moment due to sideslip, Np . The extent of 
the nonlinearity was reflected in a polynomial of degree 
sufficient to match flight and wind-tunnel data which indicated 
a nonlinear variation of aerodynamic moments and forces as a 
function of sideslip angle. A fifth order polynomial (in odd 
powers only) was proposed initially; however, as Ross pointed 
out, a cubic nonlinearity was sufficient to yield good agree- 
ment with experimental results of wing rock. Ross concedes 
that although nonlinearities existed in the dihedral effect 
and sideforce characteristic, these were considered relatively 
insignificant and thus ignored. 

An analysis follows which resulted from an intention 
to simplify the mathematical modeling proposed by Ross and 
provide a better analytical tool for use in analyzing lateral- 
directional motion at angles of attack near stall. As the 
reader will see, the analysis and conclusions presented here 
do not necessarily agree with those of Ref. [3]. 

2 . Theory 

As a means of simplification, the sideforce equation 
and the accompanying spiral divergence mode were, as in the 
hysteresis evaluation, assumed to be insignificant and thus 
eliminated from the equation set describing lateral-directional 
motion. This is not unreasonable, as flight tests involving 
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aircraft exhibiting wing rock have shown that the sideforce 
is sufficiently small as to be non-discernible to the pilot 
and has virtually no effect on the heading or ground track. 
Therefore, the equations of motion applicable to this analysis 
are the same as in the section on hysteresis, i.e., equations 



» • 

0 ' ~/Vf JS + yf - 

f lf/3 - i^js ~l f f 



Ross concluded that the term which contributed most 
to the limit cycle in view of an unstable Dutch-Roll oscilla- 
tion was a cubic nonlinearity in yawing moment due to sideslip. 
For this analysis the cubic nonlinearity was introduced as: 



Yawing moment due to sideslip = /? = /?V 

where 6 is a scalar of arbitrary magnitude which facilitated 
varying the desired level of nonlinearity for a given analysis. 
It was necessary to reduce the set of simultaneous differential 
equations in (05) to a single third order equation as a func- 
tion of a single independent variable , /& . From equation (03a) 
the cubic nonlinearity yields: 



/? “ (03c) 

Taking the derivative of this equation with respect to time 
yields : 

? - - 4 ^ 4 - 3 + *U fi - yi 

( 12 ) 
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Solving for in (12) as a precursor to elimination of the 
variable ^ 7 from the equations set yields: 

7 ? - Tty /"/ (/+3€j3')/3 (13) 



Now solve for from equation (03a) and substitute the result 
in equation ( 03 b) . 






Equations (13) and (14) can now be equated to complete the 
elimination of ^and yield a single third order differential 
equation in the single variableyf, from which the analysis 
proceeds . 



Clearing through by a minus ^ and collecting terms yields: 



+1^/? +[ / Vp(l + 3i(? x ‘) ~Nfl A + )]/f*Q 

It is convenient to introduce the following notation for ease 
of writing: 



A * -CN r * L p ) 
B - bi|S * b 3 f 
C - =1(3 * c 3 ^ 3 



where: b, = N** + N L - N L 

1 P r p p r 



b, - i N* 



C 1 ' V/S ' L p N /? 

c 3 - -«“# L p 



Thus, the resulting third order equation can be expressed in 
a more general form as : 



c it 



+ ft 






dB(g), 
</* + 






(15) 
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This is a nonlinear differential equation for which there is 



no analytical solution. Necessary assumptions were made in 
order to gain an approximate solution as discussed in the 
following . 

3 . Analysis 



motion given by (15) it was possible to proceed with an analysis 
of the motion similar to that presented by Ross, differing only 
in the degree of simplicity. Recall that the simplicity was 
gained by making the assumption that ¥ = -ft or that the air- 
craft remains on a steady track (without heading deviation) 
while engaged in wing rock. Ross, however, did not make this 
assumption but rather included the side force equation result- 
ing in a fourth order system. It was further assumed, as did 
Ross, that the oscillatory motion regardless of damping could 
be represented by the following: 



This simply says that wing rock motion can be described as a 
sinusoidal oscillation with sideslip amplitude considered as 
a slowly varying function of damping and time. The shorthand 
notation of sigma (<r*) and -fch^ta. ( d ) was used for writing 
convenience . 

In order to gain an approximate solution to equation 
(15), the Krylof and Bogoliubof (KB) method. Ref. [7] was used. 



With the resulting general form of the equation of 



/ - 0? cos u> t = COS & 



where 




(16) 
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This method has been used extensively to analyze second order 
systems with ’’weak" nonlinearities. With the exception of 
this analysis, as well as that done by Ross on the fourth 
order system, no other reference to use on higher than second 
order systems could be found. Nevertheless, there was no 
apparent mathematical restriction to applying the KB method 
to the third order system of equation (15). Moreover, the 
KB method is applicable given the assumption that the ampli- 
tudes of a vibratory system vary only slightly over a single 
cycle. Thus, the terminology often used is the assumption of 
slowly varying parameters, which is equally valid when con- 
sidering a slowly divergent Dutch-Roll oscillation. 



solution in equation (16) for the dependent variable,^, and 
derive the first, second and third derivatives, which were: 



It was necessary to take the assumed form of the 



; , - 7~ [X cos & - Co sin 6 ] 



(17a) 



= <T [ (/.*- w l ) cos e - 2 wX sin 9 ] 



(17b) 




Substituting equations (17) into equation (15) yields: 
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Equation (18) is of the form: 

= r F (Xj LU, T, COJ SJH &) 

The governing principle of the KB method says that since the 
amplitudes are considered essentially constant over a given 
cycle Cslowly varying parameters) , then the average value 
over the cycle is a sufficiently close approximation to the 
actual maximum value. This averaging process, when applied 
to equation (13), yields the following: 



£t J - 3 «V <7 - 

[^-JCoX'] = 

[X ? ~ 30J X /] = 

fu) 3 - = 

Arriving at equations 





ITT 

p S//V 6 <f 6 



if (-Trfaa l -L>')+l, t X + ej - * tr[rtjX+<j] r 1 

/ r (19a) 

77 1 X, coJ+ ifirl> J Lor x j ( 19 b) 

(19) necessarily involved the orthogon- 



ality principles of integration as shown below: 




Note that equations (19) yield the characteristic equation 
for the stable system, when the nonlinearities vanish which 
are contained in the odd subscripted variables, b^ and c^: 

i uj) - D S - Su> */] " I [^ 3 - 3U>J. Z ] 

(X h l u>) = “/? '(/. + i uj) - ( /{ f / us ) - C ; 
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which is equivalent to: 



5 3 + 6 'J% b,5 -h C , -o 

Assuming that a limit cycle occurs, the following conditions 
must hold: 

/s = 7} cos 6 

X 'b 

<rj = d 

Thus / j = 0 = damping 

17^ = D = constant 

where the subscript, s, is used to denote steady-state condi- 
tions. Hence, from equations (19), with the limit cycle 
conditions of zero damping imposed, A 5 = 0 , there is 

O = dCO; - C y - %. c 3 <77 d 

S -?/ 2. 

= 4 co s + ^ t> ? Co <Tj 

It can be reasonably assumed that a limit cycle will not occur 
for a frequency ^ = 0. The phrase "reasonably assumed" is used 
advisedly. Although it seems intuitive that such an assumption 
would be valid, some researchers, notably Parkinson [Ref. 8] 
detailed a study of the transverse galloping of a long prism 
of square section in a steady flow field in which the subject 
body possessed an unstable focus at rest. Whether this means 
that cO= 0 is a valid solution for a wing rock limit cycle, 
and thus cannot be discarded, remains largely unanswered. 
Nevertheless, such an assumption was made, as did Ross, enabling 
elimination of the U> s term between the two above equations 
yielding the final equations describing the limit cycle. 
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( 20 ) 




The frequency, U* s , for the limit cycle can then be obtained 



Recall that one of the preconditions for the limit cycle as 
proposed by Ross was that the aircraft must possess an unstable 
Dutch-Roll mode of slowly varying parameters, which is con- 
strained in amplitude by a positive cubic nonlinearity in 
yawing moment (yawing moment increasing with sideslip angle) . 
Yet there was nothing explicitly incorporated into the mathe- 
matical model, which yielded equations (20) and (21), that 
would account for this condition. Instead, it was assumed 
that a limit cycle existed due to a constrained motion in- 
stability. The modeling proceeded on this basis with the 
motion constraint explicitly incorporated by (/-h • 

This should not be interpreted to mean that a limit cycle 
will occur only if there is a divergence of the vibrating 
system initially (refer to section III. A); however, when 
applied to aircraft dynamics this is probably the most common 
way limit cycle oscillations occur, as was experienced by 
Ross and has been experienced by the author as well. 



as : 




( 21 ) 
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The effect of inducing a limit cycle via a positive 
nonlinearity on Np at high angle of attack is intuitively 
appealing. This would amount to increasing directional sta- 
bility with, sideslip, or the spring in the spring mass damper 
gets stiffer with displacement from equilibrium, until there 
is an equilibrium established between the energy contained in 
the unstable Dutch-Roll and that absorbed by the spring. 
Physically, this could correspond to an aircraft exhibiting 
a decrease in directional stability at low sideslip angles 
and high angle of attack due to a turbulent, low energy flow 
field emanating from the wing and engulfing ehe empennage 
section, thus reducing the restoring moments generated by 
stabilizing and/or control surfaces. This is often the case 
for the tactical high performance aircraft. However, as the 
tail is yawed out of the low energy flow and into a lower- 
turbulence- level , higher energy flow (higher kinetic energy 
or dynamic pressure), the directional stability will increase. 
Thus, the premise of an unstable motion constrained by a non- 
linearity in directional stability would appear to be valid. 

The F-94 airplane as originally used had a stable 
Dutch-Roll; therefore, it was necessary to destabilize the 
airplane artificially for the sake of conducting wing rock 
studies. The primary means of achieving an unstable Dutch- 
Roll oscillation consisted of altering the sign and magnitude 
of the yaw rate damping derivative, N r . The effect of this 
can be seen in figure 10, where the motion slowly diverges at 
a frequency equal to the stable Dutch-Roll. Since the 
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denominator of equation (20) is positive, and in order to 
ensure that the radical term remained positive, it was neces- 
sary to alter a second derivative such that Ab^ > c^. After 
some manipulation, the following derivatives were changed as 
indicated : 

N r = .02 
N p - -.3 

a. Predictions Based on Analytical Solution 

The following is a sample calculation of a pre- 
dicted limit cycle using equations (20) and (21) . 



A 

b 



1 



c 



1 



b 

c 



3 

3 



*1 



-(N a + L p ) = 2.4357 

N 0 + N r L p ’ N p L r = 1-7-274 

N La - N a L = 4.0916 
p p P p 

6 N p = 1.3214 6 
- €N-jL = 3.245 

P P \y 

, n / (2.4357) (1.7274) - 4.0916) 
^ - l (1 - 2.4557) (1.3214)/ 



with £ = 5 

T* s = . 1267 Rad = 7.15 Deg 

U> s = 1.3454 Rad/sec 

T = 4.67 sec. 
s 



In addition to the analog and digital computer 
programs previously mentioned, programs were likewise written 
and executed based upon the reduced third order equation (15) . 
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These are presented in Appendix D and F. These programs were 
used primarily to compare with the results based on the two 
equation system (equations 03b, 03c) to ensure consistency 
and accuracy to the maximum extent possible of the mathematical 
modeling that led to equations (20) and (21) . The same plot- 
ting program of Appendix C was used. Heavy reliance was 
placed upon the analog computer for qualitative analysis, 
where changes could be made and their effects observed imme- 
diately. 

b. Verification of Predictions 

The limit cycle predicted by the above sample 
calculation was not obtained. Rather, another limit cycle 
of lower steady state sideslip amplitude, figure 10, occurred 
almost immediately upon release from initial conditions with 
a level of nonlinearity, 6=5. 




TIME C5EO 



Figure 10. Time History of Wing Rock Limit Cycle 
and Divergent Dutch-Roll. 
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The value of the predicted <7^= .13 radians was greater than 
the initial condition /% = .08. Thus, according to theory and 
the fact that an unstable Dutch-Roll existed as shown in Fig. 
10, the motion should have slowly diverged to the predicted 
steady state values. Similarly, an initial condition greater 
than the predicted value of ^ was tried (not shown) to see 
if the motion would converge to the limit cycle; however, the 
time history was essentially identical to the first in that a 
limit cycle again developed with an amplitude very nearly the 
same as the initial conditions. Therefore, the wing rock 
limit cycle appeared to be a function of the stability deriv- 
atives, N r and Np, and the initial conditions, rather than of 
the value of the nonlinearity in yawing moment. 

The effect of the level of nonlinearity as viewed 
on the analog computer was merely to alter the frequency, with 
no change in amplitude. On the other hand, the predicted 
frequency, , agreed very well with experimental data. For 

the above sample, the predicted in terms of the period 

2TT 

(where period T = f—) was 4.67 sec, while experimental results 

s 

showed 4.7 sec. This close agreement on frequency held for 
changes in the level of nonlinearity L, as well as changes in 
the stability derivatives, as long as the requirement Ab^> c^ 
was satisfied, as previously discussed. 

Experimentation on the analog computer showed 
that the time history response behaved as if neutral stability 
existed. When considering motion without nonlinearity imposed, 
the Routh Stability Criterion can be used to determine the 
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neutral stability boundary, and it appeared at first that 
possibly the Routh discriminant had been inadvertently satis- 
fied through the manipulation of N r and N^; i.e., Ab^ - c-^ = 0, 
is directly proportional to Routh discriminant. However, this 
was not the case as Ab^ - c^ was found equal to 0.1158. 

4 . Summary of Wing Rock Induced by Nonlinear 

Yawing Moment 

Now, an obvious question arises as to the source of 
disagreement between these findings and those of Ross. The 
main difference with this analysis and that carried out by 
Ross was the simplification imposed at the beginning, i.e., 
that the sideforce equation and accompanying spiral mode were 
insignificant and need not be considered. Since the F-94 
possessed no natural instability in the Dutch-Roll at the 
angle of attack tested, the artificial manipulation of the 
stability derivatives N r and Np , as discussed above, may have 
been so arbitrary as to not correctly represent the character- 
istics previously exhibited by the Handley Page 115, and docu- 
mented by Ross. Whether these or other concepts were the 
source for the disagreement in results is unknown. Possibly 
the analysis Ross presented for the Handley Page 115 was 
peculiar to that aircraft, although intuitively one would 
tend to argue in favor of the premise that a limit cycle 
would occur for an unstable Dutch-Roll in the presence of an 
increasing yawing moment with sideslip angle. In view of 
this, it is recommended that follow-on research be conducted 
using the complete set of equations of motion, i.e., including 
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the sideforce equation and the analysis technique presented 
here to make a more conclusive determination as to the validity 
of the concept that a wing rock limit cycle can be the result 
of an unstable Dutch-Roll coupled with a cubic nonlinearity 
in N^> . Moreover, any additional research should be conducted 
on an aircraft known to exhibit wing rock motion such as the 
Handley Page 115 or preferably on an operational, tactical 
aircraft such as the Northrop T-38 or F-5. 
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IV, CONCLUSIONS AND RECOMMENDATIONS 



The existence of a lateral-directional aircraft limit- 
cycle oscillation, due to an aerodynamic hysteresis corres- 
ponding to a roll or yaw moment relay action that is dependent 
upon the sign of sideslip velocity, has been demonstrated. 

The use of state variable notation allowed the limit-cycle 
motion to be expressed in a concise and clear manner. 

Although bank angle was not a state variable in the 
analysis and hence could only be estimated from roll rate, 
the inclusion of this term by the introduction of a sideforce 
equilibrium equation is practicable for a more detailed anal- 
ysis. However, in keeping with the purpose of the analyses 
presented here, to illustrate a concept, the only equations 
used were those absolutely required to achieve that goal. 

The existence of a wing rock limit cycle oscillation due 
to a nonlinear yawing moment, as reported by Ross [Ref. 3], 
was not confirmed by the analysis presented here. The essence 
of the differences in the two analyses was in the basic assump- 
tion that the sideforce equilibrium equation need not be 
included. Such an assumption is based in large part upon 
personal flight test experience of the author, in which the 
characteristics of wing rock at high angle of attack have 
clearly indicated that yaw excursions are relatively small, 
thus adding support to the contention that sideforces be 
assumed negligible. 
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From the exhaustive experimental observations made in 
the laboratory, albeit using an aircraft that was known not 
to possess any wing rock limit cycle characteristics, serious 
doubts have been raised concerning the generality of Ross' 
analysis. While on the one hand the concept of a nonlinear 
spring, as applied to aircraft dynamics, and the "limiting" 
of the divergent Dutch-Roll oscillation, has intuitive appeal 
experimental results using analog or digital modeling consist 
ently indicated an effect only on the frequency of the vibra- 
tion, with no effect on damping. Indeed, a limit cycle was 
achieved based upon the modeling criteria used by Ross, but 
not as predicted with respect to steady state amplitude. 
However, the predicted and resulting experimental value for 
frequency of the oscillation showed excellent agreement. Of 
the limit cycle obtained, the damping appeared to be a func- 
tion of the modified stability derivatives necessary to 
artificially destabilize the Dutch-Roll of an otherwise 
stable test aircraft (F-94) rather than of the level of the 
cubic nonlinearity. 

The source of disagreement between the results obtained 
here and those of Ross remained unknown at the termination 
of this research. Sufficient time was unavailable to pursue 
the question further; however, it is recommended that follow- 
on research be conducted, when time permits, to make a more 
detailed analysis, possibly including the sideforce equa- 
tion, in an effort to determine the source of disagreement 
between the two analyses, which has led to quite different 
aircraft responses . 
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Finally, it is hoped that these analyses will make clearer 
that interpreting an actual aircraft time history as represent- 
ing a Dutch-Roll mode with exactly zero damping may be mis- 
leading, since a motion with very similar appearance could be 
a limit cycle due to aerodynamic relay, or hysteresis actions, 
or nonlinear yawing and rolling moments as posed by Ross. The 
direct problem of determining the time response of an aircraft 
due to forcing functions, be they linear or nonlinear, is 
tractable. The inverse problem of identifying the plant from 
the aircraft time histories (known as aircraft parameter iden- 
tification) is much more difficult, because of the uncertainties 
in estimating types of plant nonlinearities. 
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ANALOG COMPUTER PROGRAM FOR 
6 = -N g (l + e6“) B + N r S - N p p + AN*SIGN(§) 
P = L g 6 - L r B + L p p + AL * SIGN(B) 
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